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1. Introduction 



Although the nuclear three-nucleon interaction (TNI) is very weak in com- 
parison to the two-nucleon interaction , it nevertheless has been found that 
the binding energies of the bound three- and four-nucleon systems cannot be 
understood without taking into account the attraction caused by the TNI [0]. 
The main component of the nuclear TNI is that associated with two-pion ex- 
change, which arises from, for the major part, pion rescattering through an 
intermediate virtual A 33 resonance. The standard model for the two-meson 
exchange component of the TNI is the so-called Tucson-Melbourne model, 
which also includes p-meson exchange in addition to pion exchange [[3], |J . It 
has been found that when realistic models for the nucleon-nucleon interaction 
are employed, the Tucson-Melbourne 7r-exchange model for the TNI leads to 
an overbinding of a few hundred keV in the trinucleons f| and of 2-4 MeV 
in the case of the alpha particle ||. To compensate for this overbinding an 
additional repulsive spin-independent phenomenological TNI of short range 
has been proposed J7J. On the other hand, the Tucson- Melbourne model 
(with both 7r- and p- exchange) has recently been shown to give the correct 
binding of the tri-nucleon with the Tucson-Melbourne meson-baryon-baryon 
vertex functions ("form factors") |J. These form factors, however, differ in 
their short range behavior from those of the (older) realistic models of the 
nucleon-nucleon interaction. While this form factor discrepancy || is further 
studied in the case of the NN interaction [lOfl, there remains a need for fur- 



ther investigation of the short range aspects of the TNI. 

There is no known dynamical mechanism that would lead to a spin- 
independent short range TNI. The short range three-nucleon interaction that 
arises from scalar and vector meson exchanges with an intermediate nucleon- 



antinucleon pair were considered in ref.[ll], but were found to be both in- 
significantly small and spin dependent. We here consider another set of re- 
lated three-nucleon interactions - those that arise from pion and "effective" 
scalar and vector meson exchanges and that involve excitation of intermedi- 
ate nucleon-antinucleon pairs and iV(1440) resonances (Fig. 1). The presence 
and form of the former (Fig. la) - and more important one of these interac- 
tions - is implied by the pion and scalar and vector meson exchange compo- 
nents of the nucleon-nucleon interaction, and can be derived directly from a 
given complete model for the two-nucleon interaction. The derivation is es- 
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pecially straightforward if the nucleon-nucleon interaction has explicit scalar 
and vector meson exchanges with associated meson-baryon-baryon vertex 
functions, but the former TNI can also be constructed from "effective" scalar 
and vector meson exchanges obtained from any complete energy-independent 
model of the nucleon-nucleon interaction [|12|]. The magnitude of the latter 
TNI (Fig. lb) is less certain because of the wide uncertainty in the meson- 
nucleon-iV(1440) ("Roper resonance") coupling strengths. Both of these sets 
of three-nucleon interactions nevertheless have a very simple form, although 
spin dependent, and when combined provide an amount of additional repul- 
sion that approximately corresponds to that required and hitherto ascribed 
to the purely phenomenological spin-independent short range TNI. 



In section 2 of this paper we derive the effective 7r-scalar and 7r-vector 
meson exchange three nucleon interactions which arise from excitation of 
intermediate nucleon-antinucleon pairs and show how the potentials which 
describe these interactions can be constructed from realistic models for the 
nucleon-nucleon interaction. It is this construction method which distin- 
guishes the TNI's of the present paper from the Brown-Green M, Tucson- 



Melbourne [[3], Brazil [13|, and the Fujita-Miyazawa [T%] TNI's (the last 
being extensively used in a program of studies with light nuclear systems |l5j ) 
The latter TNI's did not attempt such a tight connection with a nucleon- 
nucleon interaction but relied on other aspects of hadronic phenomenology 
for motivation and parameter fixing. In section 3 we derive the corresponding 
three nucleon interactions that arise from excitation of iV(1440) resonances 
on the intermediate nucleon and derive the corresponding coupling constants 
from the partial decay widths. In section 4 we present numerical results for 
the contribution to the binding energy of the trinucleons, which arises from 



these interactions using oscillator and Malfliet-Tjon [16] wavefunctions as 



well as Paris and Bonn OBEPQ j8l 17 wavefunctions. An estimate that in- 



dicates that the D-state component of the trinucleon bound state contributes 
only an insignificant contribution to the net matrix element of the 7r-short 
range three-nucleon interaction is presented in section 5. Finally section 6 
contains a concluding discussion. The partial wave decomposition of the TNI 
is described in the Appendix. 
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2. Pion-scalar and -vector meson exchange three-nucleon interactions 



Realistic models for the nucleon-nucleon interaction as eg. the Bonn [fij 



Nijmegen ||19|| , and Paris |2(J models are based on meson exchange models, 
in which the Lorentz invariant on-shell nucleon-nucleon scattering amplitude 
is constructed from phenomenological meson-nucleon Lagrangians, and then 
used as an off-shell kernel in a wave-equation. In practice a nonrelativistic 
"adiabatic" approximation is also involved pi] , 22] . The important point is 
that a given model for the on-shell amplitude that is extrapolated off shell in 
this way automatically implies a model for the NN — > NNNN amplitude 
as well. This amplitude is the central component in the three-nucleon in- 
teractions that arise from excitation of virtual nucleon-antinucleon pairs on 
the intermediate nucleon (Fig. la). Thus any meson exchange model for the 
nucleon-nucleon interaction will by construction imply the presence of three- 
nucleon interactions of this type, which formally can be derived directly from 
the NN interaction model without any need for further assumptions. 



We shall here consider the 7r-scalar and 7r-vector meson exchange TNI's 
of this form. An important part of the motivation for this is the recent ob- 
servation that the off shell tt + iV — > iV+scalar and tt + iV — > iV+vector 
amplitudes are "almost" observable, in that they successfully describe most 
of the cross section for the reaction pp — > ppir near threshold [J£J. Thus 
the theoretical model for these amplitudes can be viewed as having a good 
empirical foundation. In the derivation of these three-nucleon amplitudes we 
shall use the usual phenomenological scalar ("a", "a ") and vector meson 
( "cu" , "p" ) Lagrangians, but in the end will replace the pure boson exchange 
interactions by the corresponding "effective" scalar and vector meson ex- 
change interactions that form the most important short range components 
of the nucleon-nucleon interaction. 



To construct the 7r-scalar ( "cr" )-meson three-nucleon interaction that cor- 
responds to the Feynman diagram in Fig. la we employ the ttNN and aNN 
couplings 

£*nn = ■ 7*0 > {2.1a) 

77V 
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C-aNN = g^^ijj. (2.16) 

Here 0^ is the isovector pion and <\) a the isoscalar scalar meson field and f n NN 
and g a the corresponding coupling constants. 

The three-nucleon interaction that corresponds to the Feynman diagram 
in Fig. la is obtained by retaining only the negative energy part of the 
fermion propagator for the intermediate nucleon and adding the term with 
the pion and scalar meson couplings in reversed order. In the actual con- 
struction of the operator we exploit the fact that the nucleons are nearly 
on shell, and use the Dirac equation to simplify the algebra. In this way a 
contact term operator arises, which has to be retained along with the pair 
term in the three- nucleon interaction operator. The resulting three-nucleon 
interaction has the simple form 



Here the direction of the meson momenta are taken so that they point away 
from the intermediate nucleon, that is indicated by the superscript 2 on 
the spin and isospin operators. The symbol "permutations" stands for first 
adding a term, in which the nucleon coordinates 1 and 2 are exchanged and 
then taking into account the additional (4) terms in which the intermediate 
nucleon-nucleon pair is excited on nucleon 1 and 3 in turn. 



In order to make this expression for the pion-scalar TNI consistent with a 
realistic model for the nucleon-nucleon interaction it is natural to replace the 
term — g 2 / (k 2 + m 2 a ) in the expression (2.2) by the corresponding general "ef- 
fective" isospin independent scalar component Vg(k a ) of the nucleon-nucleon 
interaction, which may be constructed using the method of ref. [P2| . In a 



similar way it is natural to replace the term (fnNN/nhr) + m V) by the 
effective isospin dependent pseudoscalar exchange potential Vp(kn) /4m% of 
the nucleon-nucleon interaction. In this way the short range modifications - 
i.e. form factors - of the simple meson exchange interactions are determined 
by the corresponding components of the nucleon-nucleon interaction model. 
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Specifically, the method of ref. ]E2J rewrites a nonrelativistic nucleon- 
nucleon potential model on the energy shell in terms of five nonrelativistic 
spin amplitudes which can be viewed as nonrelativistic limits of five rela- 
tivistic Fermi invariants. For the construction of the TNI's in this paper 
we need the isospin independent ( + ) and isospin dependent (~) Fermi in- 
variants: scalar (S), pseudoscalar (P), and vector (V). The Fermi invariant 
potential coefficients vf(k), j = A,P,V are obtained as linear combinations 
of the nonrelativistic components of a given potential. The Fermi invariant 
potential coefficients vf are functions of k 2 only which means that the un- 
derlying interactions have no energy dependence. The procedure and results 
for carrying out this program for a potential (such as Paris), which has a 



short range behavior determined in coordinate space, are displayed in [[12 
If the potential is already expressed in terms of relativistic invariants corre- 
sponding to the exchange of scalar, vector, and pseudoscalar bosons (such as 
the Bonn or Nijmegen potentials) the procedure amounts to the replacement 
(for example) 



(fnNN/mTr) 2 v P (k n ) ( f n NN /m n ) 2 A 2 - m\ 



\2 



mi 



l + kl Am 2 N kl + ml A 2 + kl ' 



(2.2a) 



where the meson-nucleon-nucleon vertex (A 2 — m 2 ) / (A 2 + /c 2 ) , is of the form 
chosen for the Bonn OBEPQ potential used in our numerical investigations. 
In either case, the short range behavior of the TNI so constructed is fully 
determined by the short range behavior of the corresponding nucleon-nucleon 
interaction, which is the issue at hand. 



In the case of an isospin-1 scalar meson exchange (ao-channel) the TNI 
that corresponds to the expression (2.2) is 



m N m^ (k 2 . + ml)(kl + ml) 

{a 2 ■ k a r l ■ t 3 + 2ia 2 ■ P 2 t x ■ f 2 x r 3 }. (2.3) 

Here g a is the aoNN coupling constant and m a and k a the mass and momen- 
tum of the exchanged ao respectively. In view of the smallness of the nucleon 
momenta in the bound states we shall not consider the non-local term in 
this TNI, which contains the momentum P 2 = (f>2 + p' 2 ) of the intermediate 
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nucleoli. In order to make the nao TNI (2.3) consistent with the models 
for the nucleon-nucleon interaction we shall replace the simple ao exchange 
interaction in (2.3) by the corresponding isospin dependent scalar exchange 
component of the nucleon-nucleon interaction: g 2 /{m 2 a + k%) — > —Vg(k a ). 

To construct the iru TNI we employ the uNN Lagrangian 

£ = ig^j^uj^ip, (2.4) 

where g u is the uNN coupling constant in addition to (2.1a). The resulting 
TNI potential is 



Vnu = - — {^^-) 2 jjy- i — Y^tttT^a^ 1 ' ^ + (permutations). (2.5) 

This interaction, which arises from the charge component of the cu-field is 
similar in form to the ira TNI (2.2), but has the opposite sign. As in the 
case of the two-nucleon interaction there will therefore be a strong partial 
cancellation between the ira and iru three-nucleon interactions, such that 
the ira typically is the stronger interaction, because of the somewhat longer 
(intermediate) range of the effective scalar meson exchange interaction. 

The irp exchange TNI has one component that arises from the charge 
component and one that arises from the spatial component. The first one of 
these has the form 



9p , fnNN ^2 ■ k w a 2 ■ k p 

m N m n (kl + ml)(k^ + m 2 pJ 



K = ~ i^^f lU ■ ^7U 7^ ■ r\ (2.6) 



and the expression for the latter is 

Vl = -1 



9p / JtvNN \ 2 ° ' fc vr 



m N m w (kl + ml)(kp + m^ 
a 2 ■ [2P 3 + ia 3 x k^f 1 ■ t 2 x r 3 . (2.7) 
The local part of the irp exchange three nucleon interaction V® is referred 



to as the "717? Kroll-Rudermann interaction" M or the "seagull" [13| in the 
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literature, and has been considered in the trinucleon before || [24], |25[. In 
the treatment of the Tucson- Melbourne TNI (an expansion of the pN — > ttN 
amplitude) V^ p is small because of a near cancellation between the inter- 
mediate nucleon-antinucleon part and the p analog Fubini-Furlan-Rossetti 
contribution to pion photoproduction. The isoscalar corresponds to the 
remaining lead term in eqs. 2.13b and 2.14a of [4a]. It was estimated in nu- 
clear matter in [4a], derived and then neglected altogether in the Brazil TNI 



D| , and is evaluated in the trinucleon for the first time in the present study. 
To make the present expressions for the iru and np exchange three-nucleon 
interactions consistent with the nucleon-nucleon interaction model we shall 
replace the bare vector meson interactions g 2 /{k 2 + m^) and g 2 p /{k 2 p + m 2 p ) 

with the corresponding isospin independent (vy(kj)) and isospin dependent 
Vy(k p )) vector exchange components of the nucleon interaction as suggested 
inref. Pfl. 



3. Intermediate N(1440) resonances 

Excitation of virtual iV(1440) (Roper) resonances on the intermediate nu- 
cleon also contributes a weak but still significant 7r-scalar and 7T- vector meson 
exchange three nucleon interaction (Fig. lb). Naturally a contribution to 
the 7T7T three-nucleon interaction also arises from intermediate iV(1440) exci- 
tation, but this is effectively included in the Tucson-Melbourne TNI, which is 
based on an off shell extrapolation of the complete nN scattering amplitude. 

To construct the na and ttuj three nucleon interactions that are associ- 
ated with intermediate iV(1440) resonance excitation we employ the effective 
Lagrangians 

f* 

CnNN* = i—^lbj^d^ ■ f(j) + h.c, (3.1a) 

£aNN* = g^*4>a^P + h.C, (3.16) 

CwNN* = Wli)*!^^ + h.c. (3.1c) 
Here tp* denotes the Roper resonance spinor field and /*, g* and g* are the 
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nNN*, aNN* and ujNN* coupling strengths, N* being an abbreviation for 
the A/" (1440). The expressions for the no and tvuj TNI potentials that are 
associated with intermediate iV(1440) excitation can then be derived in a 
straightforward way, the results being 



2&rfl* UnnI* -ho 2 - K 

m* —m N ml (k 2 + ml)(kl + ml) 



Ka = - Jl*™! J7Ti l"<? J7T S*\n.9. IT-m ^ 1 ' ^ + (permutations), 

(3.2a) 



™ m *-m N ml {kl + ml){kl + mlY h K ' ' 

Here m* is the mass of the iV(1440) resonance. It is worth noting that the 

— * 

terms that depend on the pion momentum k w in these three-nucleon interac- 
tions have the same form as the one-pion exchange interaction between two 
nucleons. 



The nNN* coupling constant /* may be calculated from the N* — > Nir 
partial decay width as 

f£_ = ^ (Ar ^ N ^ (3 3) 

Arc (m* + m N )p(E N - m N ) 

Here p is the nucleon momentum and the nucleon energy in the rest frame 
of the decaying N* . With the Ntt branching ratio of the total decay width 
350 MeV being 60% we obtain /*/47r ~ 0.031, which is somewhat less than 
one half of the corresponding value f% NN /4ir ~ 0.08 for the nNN coupling 
strength. 

The determination of the aNN* and uNN* coupling strength is associ- 
ated with considerably larger uncertainties. To obtain an estimate for the 
aNN* coupling constant we calculate it from the decay width for N* — > 

^Ms-wave aS 

777* 

Here we assume that all of the 1 = S-wave part of the tttt continuum can 
be interpreted as a broad effective a-meson. The branching ratio for this 
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decay channel is 5-15%. Assuming m a = 410 MeV at the midpoint between 
the 7T7T threshold and kinematical phase space cutoff we obtain 



*2 

^2-~0.1 (3.5) 

47T 

for T(iV* -> A^(vrvr)^= wave ) = 35 MeV. This value for gf/An is expected to 
have an uncertainty of about a factor 2. 

As the N* cannot decay into a Nu state, the coupling constant g u can- 
not be determined directly from empirical data. We shall here assume that 
9u/9u — 9a I 9<j as suggested by the constituent quark model. In the Bonn 
boson exchange model OBEPQ for the nucleon-nucleon interaction this ratio 



is 1.55 |18fl . This would then suggest that 

*2 

^- = 0.24, (3.6) 
An 

a value with which a substantial uncertainty margin also has to be associated. 

The bare meson exchange potentials in the TNI's (3.2a) and (3.2b) will 
be modified at high values of momentum transfer by shorter range dynam- 
ics in the same way as the NN interaction. To describe this short range 
modification in a way that is consistent with that of the NN interaction we 
shall introduce the same vertex factors as in the Bonn boson exchange model 
OBEPQ for the NN interaction [H| by means of the substitutions 



1 


1 




ml + kl 


kl + ml 


V ^l + kl 


1 


1 




m 2 + kl 


kl + rrii 


V + kl 

a 1 tj 


1 


1 




ml + kl 


kl + ml 


V K + kl 



(3.7a) 
(3.76) 
(3.86) 

For the form factor mass scale parameters we use the values A n = 1.3 GeV/c 2 , 
A a = A w = 2.0 GeV/c 2 . 
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4. Numerical estimates 



We shall here estimate the contribution of the 7r-scalar and 7r-vector ex- 
change TNI to the binding energy of the trinucleons with a pure 5-state 
wavefunction model. In the first estimates presented, the orbital part of 
the wavefunction is described by the harmonic oscillator and three-channel 



Malfliet-Tjon I-III model wavefunction ||16|| . The use of the harmonic os- 
cillator model is motivated by the fact that the resulting matrix elements 
can be reduced to quadrature of very simple expressions, which allows the 
qualitative features to be illuminated. The numerical values for the matrix 
elements of the three-nucleon interactions considered here that are obtained 
with the harmonic oscillator and Malfliet-Tjon wave function models turn 
out to be very similar. 

In the case of a wavefunction with only a completely symmetric S"-state 
component the radial matrix element may be expressed as 

r d 3 r r d 3 v , 
<V*>= 7^9(r,v)4v 3 (r,^ . (4.1) 



(2tt) 3 7 (2tt) 3 - 

Here 0o is the totally antisymmetric spin-isospin vector and f and v are 
differences of nucleon Jacobi coordinates, which are related to the meson 
momenta in the three-nucleon interactions as 

k n ={r+^), k a ^=v. (4.2) 

In (4.1) the function g(r,v) is the Fourier transform of the nucleon density 
function in coordinate space, which in the case of the harmonic oscillator 
model takes the form 

g(r,v)=e- T2 / 2a2 e- 3v2 / 8a \ (4.3) 
a being the oscillator parameter ^mu Q for which we use the value 0.60 fm _1 . 

The oscillator model leads to the following expression for the matrix ele- 
ment of the 7T<7 three-nucleon interaction (2.2) in 3 if and 3 He: 

< v m >= ^-A(—t r d ^- v2/2a2 -i^i 

27r 4 mjv rn^ Jo v^+mi. 
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f°° __fc 2 /2Q 2 ^ /tT« 2 T , kv . ,. 

Jo dke V feT^ ( S? ) ' (4 ' 4) 

where J 3 / 2 is a modified Bessel function. The corresponding expression for 
the matrix element of the na three-nucleon interaction (2.3) is obtained by 
the substitutions g c — > g a and m a — > m a and an overall change of sign in 
the expression (4.4), if the terms proportional to the total momentum of the 
intermediate nucleon in (2.3) are neglected. The sign change arises from dif- 
ferent isospin dependence in (2.2) and (2.3). 

The similarity in form between the 7r<7 (2.2) and nu (2.5) three-nucleon 
interactions makes it obvious that the expression for the matrix element of 
V nui is given by an expression of the form (4.4) with g a , m a replaced by g u and 
m U) and with an overall minus sign. The expression for the matrix element of 
the np three-nucleon interaction V^ p (2.6) is obtained by the corresponding 
substitution in (4.4) of g a — > g u , m a — > m p , but without any change of the 
overall sign. 

With the oscillator model density function (4.3) the matrix element of 
the N (1A40) intermediate state na three-nucleon interaction (3.2a) takes the 
form 

< VL >= L t i"f* r dve -^s 

7r 4 m* — m N Jo v A + 

The corresponding expression for the matrix element of the 7ru exchange 
three-nucleon interaction V* u (3.2b) can be obtained simply from (4.5) by 
means of the substitutions g* — > g* and m a — > m w and an accompanying 
overall sign change. 

The harmonic oscillator model of a three-body bound state has special 
properties which allow a practical implementation of eq. (4.1) as shown 
above. Other techniques must be used for more realistic wave functions 
because the triton wavefunction depends on the initial and final Jacobi vari- 
ables of the three nucleons. In general, the momentum transfer differences 
of eq. (4.2) must be supplemented by sums of Jacobi variables which, for 
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the harmonic oscillator, cancel the norm integral and therefore need not be 
considered in eq. (4.1). Further numerical calculations were made according 
to the methods described in and [3b]. The 5-state wavefunction that 
corresponds to the Malfliet-Tjon I-III interaction is given on a mesh over the 



coordinate space nucleon Jacobi variables described in |lq| . After a Fourier 
transformation of the potentials of sections 2 and 3, the three dimensional in- 
tegrals over two vector variables are straightforward to carry out numerically. 
On the other hand, in order to employ contemporary three-body wavefunc- 
tions from momentum space Faddeev calculations it is necessary to make a 
partial wave decomposition of V3 as was first done in [3b] . The final result of 
the partial wave decomposition of the potentials developed here is presented 
in Appendix A. 



In Table I we give the numerical values for the matrix elements of the 
three-nucleon 7r-short range three-nucleon interactions that are associated 
with excitation of intermediate nucleon-antinucleon pairs. The numerical 
values have been given both for the case of the oscillator model wavefunction 
and the S- state wavefunction that corresponds to the Malfliet-Tjon I-III in- 



teraction [16|. The numbers in the table correspond to the case when the 
Bonn boson exchange potential model OBEPQ (Table V, [fTgfl ) parameter- 
ization has been used to construct the "effective" pion (pseudoscalar) and 
scalar and vector exchange potentials as described in section 2 above. 



In Table II we give also for this parameterization of the TNI the nu- 
merical values for the Bonn OBEPQ and Paris wavefunction models — see 
Appendix for the technicalities involved — corresponding to the the tt-uj and 
7r-(T exchange force only (according to Table I they dominate the overall tt- 
scalar and 7T- vector exchange TNI). Table II confirms the trend of repulsion 
obtained with the more schematic wavefunction models used in the calcula- 
tions of Table I. It also stresses once more the traditional extreme results of 
the OBEPQ potential for the triton binding energy. 

All wavefunction models indicate that the S"-state matrix elements of the 
7r-short range three-nucleon interactions that are associated with intermedi- 
ate NN pair excitation are about 100 keV-200 keV (99 keV, 126 keV, 134 
keV, 194 keV respectively for the H.O., MT, Paris and Bonn model wave- 
functions), when these three nucleon interactions are constructed so as to 
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be consistent with the Bonn OBEPQ model for the NN interaction. This 
100 keV has a substantial theoretical uncertainty margin that is due to the 
remaining uncertainty in the short range behaviour of the nucleon-nucleon 
interaction. To illustrate this we have also constructed these 7r-scalar and 
it- vector exchange three-nucleon interactions from the Paris model for the 
nucleon-nucleon interaction . In that potential model the isospin inde- 



pendent scalar exchange component is much weaker than in other realistic 
phenomenological potential models, and is in fact repulsive at short range 
[EH- A s a consequence that matrix element < V na > calculated with the 
oscillator wave function model is only 143 keV (or 113 keV for a even more 
consistent calculation with the Paris wavefunction) as compared to the cor- 
responding value 392 keV obtained with the Bonn OBEPQ potential. This 
reduction of < V va > makes the net matrix element of the 7r-short range 
three nucleon attractive, when constructed from the parametrized Paris po- 
tential (~ -200 keV for the oscillator wavefunction, or ~ -101 keV for the 
Paris wavefunction). 

The difference of almost 300 keV (and a change of sign) between these 
two most consistent calculations presented here (a Bonn wavefunction and 
Bonn TNI compared to the Paris wavefunction and Paris TNI) may also in 
part be due to an inconsistency in the /x parameter of the latter combina- 
tion. The value of // labels an ambiguity in the relativistic corrections ( ie 



from v /c-expansion methods) to operators involving pion exchange |26]. The 
same continuous parameter also acts as a chiral rotation to determine the 
pseudoscalar (/i = 0) or pseudovector (/i = 1) content of the pion-nucleon 
coupling in models of the nucleon-nucleon interaction or TNI's which obey 



approximate chiral symmetry | 27fl . Fixing arbitrarily this parameter (or the 
more general parameter ft which includes a specification of energy transfer 
in the non- relativistic irNN vertex means choosing a special represen- 
tation of a unitarily equivalent class of operators. Because of this unitary 
equivalence, the observables should not depend on the value of ft, i.e. on the 
choice of representation, provided that all operators and wave functions are 
chosen consistently. The derivation of V na of equation (2.2) and of V nw of 
equation (2.5) corresponds to ft = -1, as does the one-pion-exchange part of 
the Bonn OBEPQ. Thus it is completely consistent to determine the short 
range modifications of the simple meson exchange interactions of eqs. (2.2) 
and (2.5) by the corresponding components of the Bonn OBEPQ nucleon- 
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nucleoli interaction model. The same replacement procedure with the Paris 
potential is, however, more problematic, as the parameter of that poten- 
tial has been identified to be jl = PS] . This inconsistency between the 
derivation of the TNI and the choice of short range modification is a possible 
cause of the discrepancy between our "consistent" calculations. Note that 
the harmonic oscillator NN model and the Malfliet-Tjon NN model do not 
have any pion exchange component so this particular consistency question 
does not arise. Also, the other parameter v of the unitarily linked operators 
of this paper is always v = 1/2, corresponding to no meson retardation, so 
no inconsistency can be attributed to v. 

The TNI's due to an iV(1440) intermediate state have much smaller bind- 
ing energy effects than those just discussed, a fortunate result for our goal 
of building TNI's consistent with a given nucleon-nucleon interaction. The 
matrix element of the iV(1440) intermediate state ira three-nucleon interac- 
tion V* a (4.5) for the ^-state oscillator model for the bound trinucleon states 
is 200 keV. Here we have used the parameters m a = 530 MeV and gl/^n = 
8.2797 for the mass and (nucleon) coupling constants for the exchanged a 



meson as suggested by the Bonn OBEPQ potential [18], and the short range 
form factors (3.7). The values of the meson— iViV* coupling constants are 
those derived in section 3. The corresponding matrix element of the ttuj TNI 
V* u (3.2b) is -150 keV so that the net contribution of the TNI associated 
with the iV(1440) is repulsive and ~ 50 keV. This is about one half as large 
as the repulsive contribution that is due to the TNI associated with interme- 
diate NN pair excitation. When the same calculations are made with the 
substitution of the Malfliet-Tjon wave function for the oscillator wavefunc- 
tion, the individual contributions are ~ 65 keV from V* a and ~ -45 keV from 
V* w , for a total repulsion of ~ 20 keV. This total is again small compared 
to the repulsive contribution of ~ 100 keV that is due to the TNI associated 
with intermediate NN pair excitation. Finally, in our most consistent calcu- 
lation (OBEQ wavefunction and the Bonn boson exchange parametrization 
of the TNI), the total repulsion from the TNI associated with the iV(1440) 
is about one third (+0.068:+0.194 keV) of that of the intermediate NN pair 
excitation with the Bonn OBEPQ triton wavefunction. 

In the S'-state oscillator trinucleon model and the semirealistic, Malfliet- 
Tjon NN potential wave function and the more realistic Paris potential, the 
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combined contribution of all the 7r-short range TNI's considered here is thus 
about 150 keV, which corresponds to a significant fraction of the repulsive 
TNI contribution needed to explain the binding energies of the three-nucleon 
bound states. For the Bonn OBEPQ triton wavefunction this estimate is 
enhanced, as usual, to about 250 keV repulsion. To put these results in 
perspective we remind the reader that a similar evaluation of the n — p 
" Kroll-Ruderman interaction" of Ref. [4b] with the Malfliet-Tjon NN poten- 
tial wave function yields about 340 keV repulsion when the formfactors are 
chosen similar to those of this paper. 

This estimate of the repulsive contribution due to the 7r-short range three- 
nucleon interaction would a priori be expected to be altered by the D-state 
terms (or tensor correlations) in the trinucleon bound states. A suggestive 
estimate of these alterations is provided by the comparison of the "Kroll- 
Ruderman" interaction" used perturbatively with the simple Malfliet-Tjon 
NN potential wave function in Ref. [4b] and the expectation values of the 
same interaction with triton wave functions fully correlated by two- and 
three-nucleon interactions from Table IV of Ref. [8c] . One finds a reduction 
in the contribution of this nucleon-antinucleon pair term by 20% to 40% 
from the Malfliet-Tjon result when a careful calculation with the Nijmegen 
or Paris triton wavefunction is made. In the following section we attempt to 
understand these comparisons by showing (again with the aid of the harmonic 
oscillator trinucleon model) that the effect of the D-state admixture is a 
reduction of the pair term and enhancement of the isobar term so that effect 
on the net repulsive 7r-short range TNI is very small. 
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5. Estimate of the D-state contribution 



The wavefunctions of the bound trinucleons that are obtained with real- 
istic nucleon-nucleon interaction models contain D-state admixtures of the 
order 10%. As the 7r-short range three-nucleon interactions derived in sec- 
tions 2 and 3 above have a rank-2 spatial tensor component they should be 
expected to have large SD-cross term matrix elements. In order to obtain a 
numerical estimate of the relative magnitude of this SD-state contribution 



we construct a model D-state wavefunction as |29fl 

WD = N D ^ Sij^jTi • Tj(p (j>o, (5.1) 
i<j 

where <fo is the oscillator model trinucleon orbital wavefunction 



VQ = (^L)3 e -« 2 r74 e - a V/3. (52 ) 

7rv3 

Here r and p are differences of the initial and final state nucleon Jacobi 
coordinates. The corresponding S'-state wavefunction is simply ip s = ^/Ps^p , 
where Ps is the S'-state probability. If the three-body tensor terms as S12S23 
are dropped in the calculation of the D-state normalization factor Np, it has 
the value 

N ° = ^M^~ D ' (5 ' 3) 

where Pd is the D-state probability. 

As the main source of the D-state admixture is the pion exchange tensor 
interaction it is a natural approximation to evaluate the SD cross term matrix 
elements of the three-nucleon interactions with only the nucleon pair that 
involves a pion exchange having a D-state component. The relevant SD- 
cross term density function is then 

^ 2 (r, v) = -J^S 12 (r)r 1 • rV^V 3 ^ ({U) 
where the tensor operator Si 2 (t) is defined as 
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5 12 (f) = 3(5* • f)(a 2 ■ f) - a 1 ■ a 2 . (5.5) 

With this model density function the S-D-state matrix element of the na TNI 
(2.2) can be expressed in the form 

^ T/ ^ VPdPs 91 i fnNN ^ 2 

< ' 7TCF >Si3 + D5- 



37r 4 \/5 m^v 



m-n 



2 /0 „2 Z" 00 „ _ t 2 /9 „2 k 3 



e 



f 2 + m 2 7o A; 2 + m 2 



{5^/ 1/2 (— ) - 6[k 2 + -]/3/ 2 (^) + »*/5/ 2 (£k)}. (5-6) 




The expressions for the SD matrix element of the ira^, iru and irp exchange 
three nucleon interactions that are associated with intermediate NN pair 
excitation can be obtained by making the same substitutions in eq. (5.5) as 
described after eq. (4.5). 

Assuming a 10% .D-state admixture we find the numerical values for the 
S'-D-state matrix elements calculated using the expression (5.5) to be —80 
keV in the case of the na TNI and —20 keV when the contributions of the 
7T<7, 7rao, vrp and iru three nucleon interactions are combined and constructed 
from the Bonn OBEPQ model for the nucleon- nucleon interaction ji~8|| . The 



effect of the S'-D-state cross term matrix element is thus to reduce the repul- 
sive 7r-short range contribution to the TNI that is associated with NN pair 
excitation by about 20%. 

The SD state matrix element of the na exchange interaction (3.2a) that 
is due to the excitation of intermediate iV(1440) resonances can be derived 
by the same method using the S-D-state density function (5.3). The resulting 
expression is 

. T r* . _ VPdPs 9a9a fnNNfn 

< it a >SD+DS~ 



67T 4 a/5 m* - m N m 2 
2 k 4 



J v 1 + mi J 



k 2 + ml 
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-24W 3/2 (|^) + 6» 2 / 5/2 (|^)}. (5.7) 

The expression for the matrix element of the corresponding ixuj three nucleon 
interaction (3.2b) is obtained from (5.6) by the substitutions g a , g* — > g u , g* 
and m a — > m w and a change of the overall sign. The numerical value for 
< V* a >sd+ds we find to be 94 keV, which is about one half of the cor- 
responding SS matrix element. Addition of the matrix element of the ttuj 
exchange interaction < V* u > reduces this value to 21 keV. The net effect of 
the .D-state is thus to increase the total matrix element of the 7r-short range 
TNI that arises from excitation of intermediate 7V(1440) resonances by about 
40%. 

As the numerical value for the SD state matrix elements of the 7r-short 
range three nucleon interactions that are due to excitation of intermediate 
nucleon-antinucleon pairs and iV(1440) resonances are almost equal in mag- 
nitude and opposite in sign we conclude that the importance of the D-state 
admixture in the bound trinucleon states is very insignificant in this instance. 
Hence the estimates in section 4 that were obtained with pure S 1 - state wave- 
function models and which are very similar for the wavefunction models 
considered should be viewed as fairly robust. 



6. Discussion 

The present results demonstrate that the pion-scalar and pion-vector me- 
son exchange three-nucleon interactions are important on the general scale 
of three-nucleon interactions. At the level of precision attained by calcula- 
tions with the present realistic semiphenomenological nucleon-nucleon inter- 
actions, which also contain three-nucleon interactions, this TNI has to be 
included in the calculation of nuclear binding energies. The repulsive con- 
tribution of this 7r-short range exchange TNI appears able to explain most 
if not all of the repulsion hitherto ascribed to the purely phenomenologi- 
cal spin-independent TNI of short range, which was introduced to achieve 
agreement with the empirical binding energies of the few- nucleon systems [7]. 
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The numerical values of the matrix elements of the 7r-short range three- 
nucleon interactions presented here should be realistic in spite of the employ- 
ment of pure S'-state models of the wavefunction of the bound three-nucleon 
system. By considering a schematic model for the D-state component it was 
shown that due to a number of cancellations the net effect of the D-state 
component for these TNI's is small. 

The most uncertain in magnitude of the three-nucleon interactions con- 
sidered here is that associated with the excitation of intermediate iV(1440) 
resonances. The main uncertainty in this interaction is due to the unknown 
ujNN*(1440) coupling constant. Fortunately, this interaction has smaller 
effects on the triton binding energy than the interaction due to nucleon- 
antinucleon pair terms which can be directly related to realistic NN interac- 
tions in the manner shown here. 
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Appendix — Partial wave decomposition 



While the calculations that employ the semi-realistic Gaussian and Malfliet- 
Tjon wave functions provide results that are qualitatively indicative, present 
state-of-the-art three- nucleon calculations have provided several examples of 
delicate cancellations that may be missed by schematic wave function mod- 
els. This motivated a calculation based on more realistic wavef unctions, such 
as the ones obtained with the Paris and Bonn OBEPQ potentials (in refer- 
ences || These two potentials lead to quite different answers for the 
binding energy of the triton, and different probabilities for the S and D state 
components of the wavefunction. 

The expectation value of the TNI can be calculated by means of equation 
(4.1). Although this equation is general, it requires the knowledge of g(r, v). 
Since this last function is not directly supplied by any standard Faddeev 
code, it demands extra computational effort. The algorithm that is usually 
employed with realistic interaction models is to evaluate the matrix element 
of the TNI in the partial wave decomposed basis used in the wave function 
calculation. The partial wave decomposition has been presented in a fairly 
general way in ref . [3b] . We here review its main steps, keeping the notation 
close to the one introduced in that work. 

In the LS coupling scheme the notation for the partial wave decomposed 
wavefunction is 

\pqah = \pq, [(l\)L(sl/2)S]JJ z , (tl/2)TT z } 2 

where the index 2 specifies the particle that is taken to be the spectator in 
the definition of (p, q). To calculate the TNI matrix element in this basis, we 
start by separating explicitly the spin dependence of the TNI from the orbital 
one. For this purpose the TNI is decomposed in spherical components and 
consequently in spherical harmonics. Subsequently closed form quadrature 
over the angular arguments of these spherical harmonics is performed. This 
requires expansion in Legendre polynomials of the angle dependence of the 
form factors and propagators. These expansions (in three angles) are done 
numerically, through a Gaussian mesh of at least 20 points for each angle. 
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To illustrate the technique we consider the intermediate NN n — o TNI. 
For this particular case we have (leaving out the trivial isospin dependence), 



47T. 



(a 2 ■ k^a 1 ■ K) = —KK E v / 2TTT [a(l) x a(2)f Y*(K, K). (Al) 
6 k 

Here the Y^k^, k a ) function stands for the two (coupled) spherical harmon- 
ics, which depend separately on the momenta of each meson, i.e., 



Y 1 x Y 1 



E c. 



mim2m^mi 



(A2) 



mi ,m,2 



After expressing the exchanged meson momenta in terms of the Jacobi co- 
ordinates p, q (for the initial three-nucleon state) and p', q' (for the final 
three-nucleon state), 

K = ip-p') - ^(s - 9')> 



the function Y^k^/ka) can be decomposed in coupled spherical harmonics 
of simpler arguments: 

Y l \(k 7r ,k (7 ) = 3 E E (-l) n F(l,r 1 ,r 2 )F(l,8 1 ,8 2 )(h n+n 

ri+T2=l Si+S2 = l 



X 



\p — p'\ ri+Sl \q — g*'| r 2+ s 2 

k-jr kfj 



E 

tlt 2 



n r 2 1 
si s 2 1 
k ti £ 2 fci , 



' u 000 "^000 



xY t k lt2 (p-p',q-q'), 



(A3) 



where 



F(a, 6, c) 



(2a + 1)! 



^(26)! (2c)!' 

Denote by f a {k a ) the product of the propagator of meson a by the two NNa 
couplings, including the form factor function that is introduced at the ver- 
tices. The momentum k a depends on three angles: the angle between p — p' 
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and q — q', whose cosine is x±, the angle between p and p', the cosine of 
which is X2, and the angle between q and q', the cosine of which is x 3 . 

To prepare the angular integrations we can start by doing a decomposition 
in Legendre polynomials in the angular variable x\. 

fAK)fAK)=E^^9iA\p-f\Aq-q'\)Pi 1 (^) (aa) 
h / 

The two additional Legendre polynomial decompositions — in X2 and x% 
- involve the function g^Qp — p'\, \q — q'\) of Eq.(A.4) as well as powers of 
\p—p'\ and \ q—q'\, originated by contracting the function Y£ t (p — p', q — q') 
from Eq. (A. 3) with P tl from Eq. (A. 4). The final result of the 3 decompo- 
sitions is 

ip,p',q,q') = J^dxi J_ i dx 2 J_^dx 3 P h (x 1 )P h (x 2 )Pi 3 (x 3 ) 

x \p - p'\ ai ~ t3 W- qT 2 ~ u UK)UK) . 

(A5) 

Using Eq. (A.l), each term of the J2k, V^cn gives for the orbital part of the 
matrix element, 

{{l'\')L'M'\Vt\{l\)LM) = Cf M L M ,{-) k+l+X E {-) h+h+h 

xlim E i&hU & U { f I 1 f } ( A A' X 1 
t^u I 3 4 1 J [ L L' k 

x</to(p>P / ) < /toA'(?>3 / ) , 

(A6) 



°000 °000 



with 



h+h=a { h h b J 

{A.l) 
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and 



E E C - - tl C ^F(l,r 1 ,r 2 )F(l, Sl , S2 )(-) 

ri+T2=i si+s 2 =l 

r 2+S2 { n r 2 1 



'■'2 



Si s 2 1 

£l £ 2 fc 



Hl 1 l 2 l 3 a 1 t3a 2 u(p,p', Q, Q') 



(A8) 

(«i = ri + si; a 2 = r 2 + s 2 ). 

Applying the Wigner-Eckart theorem, we extract from Eq.(A.6) the or- 
bital reduced matrix element, which ultimately we reconnect to the spin one, 
generating the desired final result. 
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Table I 





H.O. 


MF I-III 




0.392 


0.457 




-0.012 


-0.010 




-0.301 


-0.337 


v c 


0.020 


0.016 


TOT 


0.099 


0.126 



Matrix elements in MeV of the 7r-scalar and 7r-vector exchange interac- 
tions which involve an intermediate nucleon-antinucleon pair. The three- 
nucleon interactions are constructed from the Bonn boson exchange model 
OBEPQ for the NN interaction The matrix elements are obtained with 
S'-state oscillator and Malfliet-Tjon I-III wavef unctions. 



Table II 





PARIS 


OBEPQ 




0.551 
-0.417 


0.836 
-0.642 


TOT 


0.134 


0.194 



Matrix elements in MeV of the n-a and tt-u exchange interactions which 
involve an intermediate nucleon-antinucleon pair. The three-nucleon inter- 
actions are constructed from the Bonn boson exchange model OBEPQ for 
the NN interaction |1| . The matrix elements are obtained with the 5-state 



components of the Paris and Bonn (OBEPQ) wavefunctions. 



Figure caption 



Fig. 1 (a) 7r-short range exchange three-nucleon interaction that in- 
volves an intermediate nucleon-antinucleon pair, (b) 7r-short range exchange 
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three-nucleon interaction that involves excitation of an intermediate iV(1440) 
resonance. The wavy lines symbolize scalar and vector meson exchange. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9503013v2 



